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PREFACE 


W HILE lecturing in 1913 on the calculation of self and mutual 
induction, the writer noticed that the integral in the case of co- 
axial circles was of such a form as to make the application of the 
arithmetico-geometrical means obvious. An extremely elegant formula 
{Proc, Boy, Soc. A, Vol. c. 1921, p. 63. Appendix, this book, No. 35) 
when put to numerical test justified the method as one well adapted to 
computation. The calculation of more complex cases of self and mutual 
induction led to several interesting formuLae for the complete third 
elliptic integral. With a view to making a complete study of the group 
of recurrence formulae associated with the a.g.m. scales, further work 
was postponed until a vacation in 1921, when many of the results in- 
cluded in this book were obtained, and were thought by the ‘writer to 
be new. On later consulting the Collected Works of Legendre, Gauss 
and Jacobi, the writer found that he had rediscovered many known 
formulae, but that quite a number, likely to be of service in computation, 
appeared to be hitherto unknown. In the summer of 1922 the writer 
decided, as a vacation task, to include in a single monograph the entire 
theory of elliptic functions associated with the a.g.m. scales, thus com- 
pleting, and to some extent adding to, the work of Gauss left unpublished 
after his death, and placing before mathematicians in accessible form 
a mode of approach to Elliptic Function Theory directly related to the 
art of machine-computation. It is interesting to note in this connection 
that this subject was to have formed the content of the third volume 
of Halphen's Fonctiom EUiptiqties, unfortunately left incomplete at the 
time of his death. 

While the present monograph contains many new formulae and 
methods of computation, no claim is made as to novelty in fundamental 
analytical treatment. References to more important books and memoirs 
consulted by the author are given in foot-notes. A complete bibliography 
is not given since the reader will find exhaustive references in the RoycU 
Society Index^ 1800-1900, Vol. I, Mathematics, and after 1900 in t^ 
mathematical volumes of the Internatioml Catalogue 'tf Sdenti^ 
Literature, 
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league, Professor James Harkness of McGill University, for reading over 
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ON THE DIRECT NUMERICAL CALCULATION 
OF ELLIPTIC FUNCTIONS AND INTEGRALS 

SECTION I 

INTRODUCTION 

Although the numerical evaluation of elliptic functions and integrals 
is dealt with in many of the standard treatises, it cannot be said that 
the use of tables is entirely satisfactory in dealing with computations 
required in several branches of physics and astronomy. Many of the 
integrals appearing in these problems can only be expressed as some- 
what complicated expressions involving the complete elliptic integrals 
of the first and second kinds, and tlie corresponding functions of com- 
plementary modulus. These formulae give rise in many cases to 
differences of nearly equal quantities : the numerical calculation of such 
expressions to a sufficient number of significant figures is thus often a 
matter of some difficulty, requiring tedious and sometimes uncertain 
interpolations. 

While the numerical calculation of the Theta-functions is extremely 
ra])id by the use of the highly convergent series involving the nome </, 
these are not usually the functions which make their appearance directly 
in physical or astronomical problems. The necessary transformations 
required to express the results in terms of the (/-series require in many 
cases somewhat conqjlicated analysis, while numerical computation by 
this method necessitates the exact calculation of a largo number of 
auxiliary quantities*. 

The forthcoming tables in course of construction under the auspices 
of the British Association t are so arranged as to allow a large number 
of functions to be evaluated, and will very greatly aid numerical work 
in many practical problems involving the use of elliptic integrals and 
functions. While sufficient for most purposes, the use of these tables 
will still involve considerable interpolation when it is required to take 
out values corresponding to arguments and moduli expressed by six or 

* On this point see Bosa, E. H., and Grover, F. N., Bulletin 169, Bureau of Stan- 
dards, WashinKton (1913), pp. 67 and 73; Nagaoka, H., Tok. ColU Sc. J. 27 (1909), 
No. 6, and Jour. Tok. Math. Phys. Soc, 4 (1908), p. 284; 6 (1911), p. 10; Nagaoka and 
Sakurai, Scientijic Papers of the Institute of Physical and Chemical Research, Vol. ii. 
Dec. 1922, pp. 1-67; Olshausen, G. R., Phys. Rev. 51 (Dec. 1910), pp. 617-636. 

t See report (by Sir A. G. Greenhill), Brit, Ass. Report, 1912 (Dundee), pp. 
39-66. 
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INTRODUCTION 


seven significant figures, as is now necessary in several modern applica- 
tions of elliptic functions. 

For several reasons it is advantageous to have available a direct 
method of computation, independent of auxiliary tables, based on the 
use of the modern high-speed calculating machine. The writer has 
therefore thought it desirable to set out in detail the method of 
calculation based on Landen^s quadric transformation and the use of 
the scale of arithmetico-geometrical means, recapitulating briefly the 
original developments of Lagrange, Legendre, Gauss, Jacobi and others, 
extending them in several directions and keeping in mind the ultimate 
use of the calculating machine for numerical work. 


SECTION II 

HISTORICAL NOTE ON LANDEN’S TRANSFORMATION AND 
THE VARIOUS SCALES OF MODULI AND AMPLITUDES 

The importance of realizing rapid and accurate methods of calculating 
the elliptic integrals now denoted by 

u = F{<f>y k)- ( k) and k)d^ •••(!)> 

yo yo 

where A (<^, k)- sin* </>), 

was first remarked by Euler*'^ (1766), although it was not until several 

years later that Ijandenf (1775) discovered in geometrical form tlie 

transformation which forms the basis of the existing methods of the 

numerical calculation of the elliptic integrals. 

A method of successive transformations for the ultimate reduction of 
the algebraic forms of these integrals to elementary integrals was pub- 
lished by Lagrange t hi 1784-85. This memoir contains an exposition of 
the scales of arithmetico-geometrical means, together with two types of 
algebraic transformations corresponding to the increasing and diminishing 
amplitudes of Landen's trigonometrical forms: it is shown that the 
given integral must lie between two limits which may be made to 
approach each other as closely as may be desired : the limits thus ob- 

^ Euler, Novi Cornm. Acad. Sc. Petrop. vol. x. 17(56. 

t Landen, Phil. Trans. Pay. Soc. vol. lxy. p. 283, 1775; MathenuUical Memoirs^ 
London, 1780. 

X Lagrange, “Sur une nouvelle m^thode de calcul integral pour les diifdrentielles 
affect^'cs d’un radical carr(5 sous lequel la variable ne passe pas le quatri^me degr6,” 
M^m. de VAcad. roy. dcs Sc. de Turin, t. ii. 1784-5; (.Envies (Gauthier- Villars, Paris, 
1868), t. II. pp. 253-312. 
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tained for the integrals of the first and second kinds are practically 
identical with the limits and series obtained about the same time by 
Legendre. 

Landen’s transformation was applied by Legendre* to the numerical 
calculation of the elliptic integrals, and furnished the method of com- 
putation by means of which the latter^s well-known tables were con- 
structed. The scales of moduli and amplitudes employed by Legendre 
are briefly described below. 

If we form successive moduli and amplitudes according to the 
recurrence formulae 

~ (1 — )/(l + A’yt ) and tSiii ~ kji tan 071 •••(2), 

Landen’s transformation in trigonometrical form leads to the result 
^ (0n+lj ^n+l) ““ + ( 3 ). 

As n increases, the successive moduli kn converge to zero, and the 
complementary moduli kn to unity, while the amplitudes increjise so 
that 07i/2" tends rapidly to a definite limit : ultimately wo may write 
kn) ~ 0w and thus obtain the result 
h) ^ <t>nl"2^y where Ar= Jw ( 1 + X*i) (1 + . . . (1 -i- kn) 

.........( 4 ). 

The result of calculating successive moduli and amplitudes in reverse 
order is equivalent to calculating in forward order these quantities 
according to the recurrence formulae 

kn+i = 2jknl(l + kn) and sin (20„+i - 0n) = sin 0^ .--(S). 

This transformation then gives 

h+l) ' J (I +h) ^n) (fl). 

As n increases the successive moduli kn^i converge to unity and the 
amplitudes to a limit we thus obtain, when n is sufficiently large, 

/ «ec 0 n#n-logtan(i 7 r+ ( 7 ), 

giving finally 

^(«Aoi h) = >J(Jc\hk^ ... kn)IJh • log tan (Jtt 4- •••(fl)* 

Legendre also developed for purposes of numerical calculation extremely 
convergent series for k) which he employed in the tabulation of 
this function. 

* Legendre, “M<^moire sur les integrations par arcs d’ellipse” and “Second 
memoire, etc.,” M€m. da VAcad. de% Science» de Panti, ann. 1786 (Pans, 1788), 
pp. 616-643 and 644-683 ; Traite des Fomtiom ElliptiqiieSf Paris, 1825, t. i. 
p. 79 at seq. 

1—2 
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The method discovered independently by Gauss and called by him 
the “algorithm of the arithraetico-geometrical mean” originated in 
connection with the evaluation of complete elliptic integrals arising 
from a problem in attractions required in planetary theory*: in a note he 
mentions that the results were obtained by him as part of a more com- 
prehensive theory, independently of the results of Lagrange and Legendre 
with which they are closely connected. Gauss employed a trigono- 
metrical transformation along the following lines. 

Commencing with two numbers />»o), successive numbers {an, bn) 
are calculated from the recurrence forrinilae, 

^^M+l ” \ {a-n, ^n)> bn-i-i — tj (anbn) 

Writing co.s“ + /V sin^ </»»), Gauss employs the recurrence 

formula 

tail <f>Ti+i = tan <f>n , 

or tan -- tan • s/(bn + ^n)l sj{(ln + A,,) . . .(lo;. 

From the relation 

tan <l>n = (Aj Ag ... A,,...)/(«o«i ... ...) . tan <^o 

we have ultimately, since the <^’s tend to a limit <t^n> 

* 

d<l><,l\--=<f>nla«, 

and hence 

«„j^‘’rf<^„/A(, = tan-‘[(A,A 2 ... A„ ...flr„...).tan^„] ...(H). 

Extremely convergent series for the elliptic integral of the second 
kind were also obtained by Gauss. The application of this method to 
the numerical calculation of the elliptic integrals was considered by 
Jacobi some years later t. 

In the course of researches carried out between 1797 and 1808, but 
not published until after his death considerably later J, Gauss developed 
the theory of elliptic functions to a remarkable extent and derived a 
large number of formulae which, at a later date, were obtained inde- 
pendently by Abel and Jacobi and which are now associated with the 

* Gauss, “ Determinatio Attractionis quam in Punctum quodvis Positionis datae 
exerct*ret Planeta, etc.,” Comm. Gott. Soc. Beg. Scient. iv. 1818; Gcs. Werke 
(Gottingen, 1866), Bd. iii. pp. 331-3.55. The note mentioned above is dated Feb. 9, 
1818 ; Ges. Werke, Bd. iii. pp. 357-360. 

t Jacobi, Fnndamenta Nova, 1829, sections 38 and 52 ; Ges, Werke (Berlin, 1881), 
Bd. I. pp. 154 and 203. Also “ Numerisohe Berechnung der Elliptischen Functionen,” 
Crelle, Bd. xxvi. pp. 93-114; Ges. Werke (1861), Bd. i. pp. 345-368. 

J Gauss, Ges. Werke (1866), Bd. iii. pp. 361-403 and 433-469. 
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theory of the Theta-functions. In one of these postliumous papers, Gausi 
employs the recurrence formula tan2<^n= \/(^«+i/»n+i) tan which 
is much more convenient for purposes of numerical calculation than (10). 
The last mentioned method is considerably simpler than that of 
Legendre's scales of moduli, especially when a calculating machine is 
available, although in theory these two methods of computation are 
practically identical. 

Although a few formulae based on the scales of arithmetico-geometrical 
means were briefly considered by Jacobi t, no attempt seems to have been 
made by him to apply the results to actual computation. 

As the method of evaluating the elliptic functions by the use of the 
A.G.M. scales is only briefly dealt with in existing text-books and does 
not seem to have been developed to the extent it deserves in various 
memoirs which have appeared on the subject in recent years, a tolerably 
full account of the subject is given in the present book, together with 
a brief outline of the main analytical steps. Several new formulae have 
been obtained by the writer which facilitate the application of this 
method of numerical computation and considerably extend its scope. 

From the point of view of pure analysis, the modern Weierstrassiaii 
notation presents many advantages in regard to elegance and symmetry, 
a superiority which no longer holds, however, when it is required to 
proceed to numerical evaluation. As the main point of the present 
paper is concerned with improvements in methods of computation, the 
older Jacobian notation has been adhered to as far as possible, especially 
as it expresses most readily the elliptic functions and integrals which 
arise from integrations involving circular and hyperbolic functions. The 
results are set out in such a manner as to provide a fairly complete 
compendium of formulae intended to be of service to the computer in 
reducing to numerical results such expressions involving elliptic functions 
and integrals as may arise in physical and astronomical problems. 

* Gauss, loc. cit. p. 388. 

t Jacobi, Ges. Werke (1881), p. 357. 

X Cayley, A., Elliptic Functions , 2nd Ed. (G. Bell and Sons, London, 1895), 
pp. 326-338 ; Enneper-Miiller, Elhptische Functionen (2nd ed. 1890), pp. 361-364 ; 
Tannery-Molk, Fonctions ElhptiqueSf t. iv. 1902, pp. 269-275. According to a foot- 
note in Halphen’s Fonctions Elliptiques, t. ii. 1888, p. 310, the subject of the 
arithmetico-geometrical mean was to have received special attention in the author’s 
third volume, unfortunately incomplete at the time of his death. 



6 


SCALE OP ABITHMETICO-GEOMETBICAL MEANS 


SECTION III 


ON THE SCALE OF ARITHMETICO-GEOMETRICAL MEANS 


In forming the scale of arithmetico-geometrical means, we start with 
the positive numbers {ao, bd), of which ao is the greater, and form 

successively 
== K«o + * 0 ), 

Oa = + ^1). 

bi — iJ{aQbd), 
b^ — \l {ai bd), 

Cl — ^ (®0 ~ ^ 0 ) 

, ...(12). 


bn — »J{an-\bn-l) 

” ^{Hfi-l ~~ bn^d) 

We obviously have the relations 





Cn - <In — 

(13). 


bn — ^n+i “ ^n+l> 

an^n — ^ J 

The a's and i’s 

tend to the same limit denoted by 



II 

lim an 

(14) 




with extraordinary rapidity, even when and are initially numbers 
of very different magnitudes*. 

We notice in passing that satisfies a homogeneity relation 

of the form 

€3f(ao, b^ = M{faoy c^o) (15)> 

€ being any number. The array of numbers (12) will be referred to as 
the “ scale of arithmetico-geometrical means (a©) or more briefly 
as the A.o.M. scale bo). Also, the limit M{ao, bo) will be denoted 
by an, as long as by so doing no ambiguity is involved. 

If we calculate in the same way an array of numbers commencing 
with {ao=ao, bo =Co- ^{ao -bo)}, we derive what will be called the 
'‘complementary a.o.m. scale {ao, symbols being in this case 

denoted by accented letters, and the limit M {ad y bd) by . 

* The usual arithmetical process of extracting square roots is easily adapted to 
most calculating machines of modern type. An account of these is given in 
Modem Instruments and Methods of Calculation^ edited by E. M. Horseburgh 
(London, G. Bell and Sons, 1914) ; also by d’Ocagne, Le Cakul Simplifi€ (Gauthier- 
Yillars, Paris, 1905). In carrying out the process the work is considerably shortened 
by keeping in mind the rule that if the first p digits out of the number n 
required in the square root have been obtained by the usual process, the next 
jp - 1 digits can be obtained by division only, with a possible error of 1 in the last digit 
(Chrystars A Ziifc 6 ra, Part i. 6 th ed. 1904, p. 210). To illustrate the rapid convergence 
of the A.G.M. scales. Gauss works out an example starting with ao=:l and : 

a^ and 65 are represented by the same number to fifteen significant figures. 
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It is obvious that the operations of the A.a.M. scale may be carried 
out backwards : if the entries are denoted by negative suffixes, they 
are connected with those of the complementary scale by the relations 
a-n = , h-r, = (16). 

These relations evidently hold good when accented and unaccented 
symbols are interchanged. As a result of these relations it will be 
seen that there is no advantage to be gained from considering any other 
than the a.g.m. scale of positive suffixes (an, ^n) and the corresponding 
complementary scale (a^ , bn\ 


SECTION IV 

LANDEN^S SCALE OP INCREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second 
kinds, modulus k. 

If we write An= v/(an^ cos^^»4- 6/ sin* the recurrence formula 

tan(<^n4.i-^„) = (6nK)tan<^n (17) 

may be written in either of the forms 

sin (2<^n “ <#>n+i) - (Cn+l/«n+l) sin <^n+i, 

cos (2<^n ^n+i) ~ A,i4i/an+l (l^)* 

From these we derive 

^n+i ^n+i cos ^n+i “ and An+1 •— Cn+1 cos <^n+l “ ^n^n/An •••(19), 
which give on differentiation, 

d<t>ol^o — h “ i d<f)2l A 2 — ... = ( J )” d4^n / An — (20). 

From the relation 

AVi + Cn+i An+i COS i>n+i=^-h^n+h(lnbn 

we obtain, making use of (20), 

^nd<bn “ ^n+id<l>n+i COS ^n+l^^n+l “ ^nbnd<l>nl^n •••(21). 
Subtracting from each side the identity 

and<hnl^n "■ ^\+ld^n+il^n+l ^ \(^nd4>nl^n ” 0>nbnd^nl^n • • ^( 22 ) 
we obtain 

(An — anl^n^d<f>n (An+i 0'\+il ^n+i^d<f>n+l 

— Cn+1 COS <^n+l^^n+l JCn^^^n/An • • •(23). 
Writing down a series of such equations, commencing with n = 0, 
and noting that as n increases, (A,^ - tends to zero in the limit. 
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we have, on integrating, 

r<t>o 7 

I I Ao = Cl sin + C 2 sin <^2 + ••• +c„sin<^n+ 

Jo Jo 

- i (c„'“ + 2c, “ + 4 c 2 * + . . . + 2”c„“ + . (24). 

If we now construct the a.g.m. scale (ao = l, = calculate 

successive amplitudes commencing with = we have Ao = A(<^, i). 
As n increases, the a*s and //s tend to the same limit while the 
angle <#>n/ 2 ’‘ tends to a finite limit. We thus obtain from ( 20 ) and (24), 


u - k) = (IK) (0J2-) (25), 

k) - k) = (ci sin + Ca sin <^2 + • • • + sin + • • •) 

- i(co^ f 2ci^ + 4 c 2- + ... + 2"Cn^ + ...)F((I>, k) (26). 

If we commence with <^o = ^tt, we have 

= TT, <^2 - 27r, <^3 == 47r ... = Jtt, 

so that K = W/anj and 


(A"- F)IK^ i (co® + 2ci2 + 4 C 2 ' + . . . + 2’‘c„^ 4 . . . .) (27). 

The series (24) and (27) converge with extreme rapidity, resembling 
in this respect the (/-series with which they are closely connected*. 


Jacobi’s integral Z (</>, k) is defined by the relation 

(<^, k) = k) - {EIK) F{^, k) (28) 

or, if the angle ^ is connected with the argument u by the relation 

sin <^ = sn (ti, k\ 

Z{ti,k)= f (!Li\^{ujk)dti--{EIK)u (29). 

Jo 

From (26) and (27) we derive the expansion 

Z(<f}, k) Cl sin <^i + C 2 sin <^0 + ••• + Cn sin </>,!+ (30). 

In terms of Jacobi’s Theta-function 0 (u, k), 

Z(u,k)du = J^^ Z'(<f>, k) d</>/ A(<f>,k) = log [S(u,k)/&(0,k)] ...(31) 

where 0 (0, k) = J(2k'Klir) (32). 


* These relations are discussed by Legendre {Traits des Fonctiom Elliptiques, 
t. III. Deuxi^me Supplement, 1828, p. Ill et seq.]. On the connection between the 
A.o.M. scale and the modern transformation theory of the Theta- functions, see 
Tannery and Molk, Fonctions Elliptiques, t. iv. 1902, note 2, pp. 269-273. A develop- 
ment of the theory of the Theta-functions based on Legendre’s scales of moduli is 
outlined by Bichelot, Correspondenz mit Herrn Professor Schroter (Kdnigsberg, 
1868). See Examples 19-24 of the Appendix. 
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From (20) we have = so that the series (30) 

enables us to write 

r Z(</», k) (^<^/A(<^, k) = i (Cnjr) sin 

/o 1 Jq 

On integrating and making use of (19) and (18), we obtain 

r^n 

/ sin 4>nd<l>JAn = (l/c») log (a„_,/A„_i) = (l/c„) log sec (2<^„-3 - 

JO 

leading to the expansions* 

rziu, k) du = i (|)» log (33), 

JO 1 

or 

log [® («, i-y® (0, ^•)] = I log (fflo/^o) + i log sec (2(f,„ - <l>j) 

+ ^ log sec (2<iii - <^) + . . . + (0"+* log sec (2<^„_, -</>„) + (34). 

In the series (34) it is more convenient to write 
«o/Ao COS (<l>i - <^o)/cOS 00 

for purposes of numerical calculation. 

(ii) Calculation of the Jacobian Functions sn (w, i), cn (w, k\ dn (Uy k) 
etc, in terms of the argument u. 

When the argument u is given and the modulus k is known, it is 
only necessary to compute the a.g.m. scale (</o=l> ^o = ^) *** 

value of n that Cn is less than the small quantity determining the order 
of accuracy of the calculations. Since u ^ F(f>^ k), equation (25) enables 
us to calculate in circular measure the angle 0n from the relation 

0„ = 2'"«hM (35). 

From the recurrence formula (18), 

sin (20n-i - 0n) "" {Cnjan) sin 0„, 

we are enabled to calculate successively the angles 0n-i ••• 0i and 
finally 0oi in terms of which 

sn(Uy k) = smd>Qy cn (w, A;) = cos 0o> di\(u, k) -=\ = cos<t>o/cos(i>i — i>o) 

(36). 

Z(uy k) may then be computed from the series (30), and ® («^, k) from 
formulae (32) and (34). 

* Formulae (30), (33) and (34) are practically the same as those given by Legendre 
(Traiti des Fonctions Elhptiques, 1. 1 . section 90), and are contained implicitly in 
formulae given in Gauss’ posthumous papers. Formula (34) is equivalent to that 
given by Jacobi {Crelle^ Bd. xxvi. pp. 93-114; Oesammelte Werke (1881), Bd. ii. 
p. 357). 
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SECTION V 

THE HYPERBOLIC SCALE OF INCREASING AMPLITUDES 

(i) Calculation of the Elliptic Integrals of the First and Second 
kinds^ modulus k\ 

Writing if^n for throughout the formulae (17) to (24), we obtain, 
instead of (17), 

tanh - <l>n) = (bn/an) tanh (37), 

which may also be written 

sinh (2<^n “ <#»»+i) = (cn+i/a«+i) sinh cosh (2<^„ - <f>n+i) = ^n+iK+i 

(38). 

If we denote A,^ = J(an coslP <l>n - bn sinlP <^„), it is easily proved, as 
before, that 

d<l>J\ = i#,/A, = = . . . = Qfrd^hnl^n = (39) 

from which it follows that as n increases, we have in the limit, 

[ " ^^<#>o/Ao = (IK) (<i^n/2^). 

JQ 

If we now write 

sinh <;(>o = tan </> (40) 

it is easily proved that 

Aq = A (</», X;')/cos = d<hl^ i^yk') (41), 

so that 

u = F(<l>> k') = (IK) (<^n/2") (42). 

Formula (24) becomes, 

7 

/ l^od<t>o - I d<l>J^o = (ci sinh + Cg sinh <^2 + • • • + sinh + • • • ) 
Jo Jo 

-K^o" + 2ri^ + 4c2^+...4.2”Cn^+...) •••(43). 

yo 

It is easily shown that 

f = f\ (</>, k') d<htcos^ ^ 

Jo Jo 

= { dn® (m, X:') dujcv? (w, A:') = f dn^ (m, ^) 
yo yo 

or in terms of the ^T-function of imaginary argument, 

{ Aod<^o = (lli)Z{iUf k) + {EjST) u, 

Jo 

Substituting in (43), and making use of (27), we derive the series, 
corresponding to that of (30), 

(lli)Z(iUy X:) = Cisinh<^i + C2sinh<^2+ ... +CnSinh</>n+ ... .*.(44). 
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The following formula, corresponding to (34), easily follows by 
integration, 

log [© (m, k)l® (0, k)] = 4 log (aj\) + 2 log sech (20n.i - <^n) 

(45) 

in which we may conveniently write aj\ = cosh (<^i — <^o)/oosh <#>o. 

From (42) we easily deduce the relations, 

(l/i) sn (iu, k) = sn (m, k')/m (Uy k') — sinh 
cn (zw, k)= 1 /cn (Uy k') = cosh </>o 

dn (iuy k) = dn (w, kycn (^^, ^') = ^ cosh <^o/cosh (<^i - <^o) 

(46). 

Making use of Jacobi’s well-known relation in the form 
Z(Uy k') = sn (m, k') dn (^^, k')lm {Uy k') - Trul(2KK') + iZ (iUy k) 

(47) 

we have from (28) and (46) 

B(4>yk^)-^(£;'ir)F(<t>yk') 

= sinh <^o/cosh (<^i ~ <t>o) - ‘?rjP(<^, k')j(2KK) + iZ (m, A;). 
Noting that Legendre’s relation gives 

ElK + F/K'-l^i^liKK') (48) 

we obtain from (27) 

E{<t>y k') = sinh </>o/cosh (<^i — </)o) + (1 — E/K) k!) -f iZ {iu, k)y 
or, finally, making use of (27), (42) and (44), 

E{<f>, k') = sinh ^o/cosh (<^i — </>o) 

+ 2c, ^ + - + 2"c„^+ ...) (IK) (<^„/2’^) 

- (ci sinh + Ca sinh <^2 + • • • + sinh </>„+•••) • • -(49). 

(ii) Calculation of the Jojcohian Functions sn (w, A:'), cn (w, A;'), 
dn (i^, A’), etc, in terms of the argument u. 

Compute the a . g . m . scale («©= 0o = k') and calculate <hn from the 
relation ^n^^'^anUy n being such that is less than the small quantity 
determining the order of accuracy of the calculations. The hyperbolic 
recurrence formula 

sinh (20n-i - ^n) = (Cn/«n) sinli </»„ 

enables us to calculate successively ... a,nd finally <^o, in 

terms of which we have from (46) 

sn {Uy K) = tanh <^o , cn (m, K) = 1/cosh <#>o» dn {Uy k') = l/cosh (</>! ~ </>o) 

(50). 

Z{Uy K) may then be calculated from (47) and (44). 
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(iii) Calculation of the Jacobian Functions of imaginary argument, 
Hn {iu, k), cn (iu, k), dn (m, /:), etc. 

The procedure outlined above, together with equations (46), obviously 
enables us to calculate the Jacobian functions of imaginary argument 
and modulus k, while formulae (44) and (45) serve for the functions 
Z{iu, k) and ®{iu, k). 


SECTION VI 

GAUSS’ SCALE OF INCREASING AMPLITUDES 
(i) Circular Scale, 

In connection with the a.o.m. scale, Gauss made use of the recurrence 


formula* 

tan Xk+i = -s/(«n+i/^«+i) tan 2xn (51). 

It is easily proved that the above formula is connected with Landen’s 
scale (17) by the relation 

tan x,i = sjihnjan) tan 4>n (52) 

or 

^71 + 1 “■ — 2Xn-i (53). 


In the limit, as n tends to infinity, (bjan) approaches unity, so that 
Xu and both tend to the same limit. It thus follows from (25) that 
if we form the a.o.m. scale (ury=^ 1, b^^k'), and starting with Xo deter- 
mined by tiinxo - V(^>,>/oo) tan sf>, compute successively Xi» X2» ••• X» ••• > 
we have 

mC-(lK)(x./2") (54). 

From (52) and (53) it is a simple matter to derive successively 
. for use in the forimdae of Section IV. When these angles 
are required, it i.s, however, simpler to make direct use of Landen’s 
recurrence formula (17) or (18). 

(ii) Hyperbolic Scale. 

It is evident that formulae (51) and (52) may be written in the 


h)q)erbolic form 

tanh Xu^^i =" tanh 2xu (55), 

tanh Xu = Jibnla„) tanh <^« (56). 


* Gauss, Oes. Werke (1866), Bd. iii. p. 388. To compare Gauss’ results with 
those of the present paper we note that 1* appears 

(p. 392) that formula (26) of the present book was also obtained by Gauss, inde- 
pendently of Legendre (foot-note on p. 9). 
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As n increases indefinitely, Xn and <f>n tend to tlie same limit, so that 
if we write 

tan <l> = sinh or sin = tank <^o (i>7) 

it follows from (42) that 

F(<l>, k') = (l/a„) (.^„/2“) = (l/a„) (x„/2’') (M). 

Noting that, as before, 

4>n+l-<t>n = "2Xn-l (5S) 

it is a simple matter to derive successively <^,, <^ 2 , ... 4 >„ ... for use in 
the formulae of Section V. In such cases, however, it is simpler to 


make direct use of the hyperbolic recurrence formulae (37) or (38). 

(iii) Calcukitum of K', E' and the nome q -e 
If we write <^= Jir, it follows from (57) that tanh 1, and from 
(56) that tanh Xo = >J(f>Ja^. From (55) we deduce tliat 
tanh2xo tanh 2xi = /'.j/ffl., ••• tanh 2xu =- 

and tluis 2x„ =■ i log {(«-„+, + ,,)/(«'„+, - 5,,*,)} "= i 1%' («..+i.>/f»+ 2 ) 

or (x»/2") -= (i)”^’ log (rt,H 2 /n,vi)- 

It follows from (58) tliat as n is indefinitely increased, (i)" log(rt„/c„) 
tends to a finite limit given by 

a„F(\v, k’) = l^K'j K i log ( 1 A/) (),)” log (^/„/e„) . . (60). 

Making use of the relation (13), rt„c„ we have 

log (rt„/c„) = 2 log (2 <t„/«,_,) + 2 log («„-,/c„-,), 
which gives by successive application, 
log(a„/c„) - 

2 log (2((„Attt.,) + 4 log (2«„_,/rt„_2) + . + 2’* log (2rt,A/„) + 2" log (a„/c„). 

Inverting the order of the series, we obtiiin 

(D” l0g(rtnA-«) - log («„/c„) + l0g2 {1 + t + J + ... + (I)”-!! 

+ log (a,/rt„) + i log («,/«,) + .. f (’)’•-' log («„/«„.,). 
As n is increased indefinitely we obtain the series 

«„A'' = -tlog(4a,/c,)-i;(J)''log(«„A/„^.) (61). 


Making use of Legendre’s relation (48) in the form 
E' ^ (I -EjK) K' + a„ 

and using (27), we obtain the following useful formulae, 

A” = \ ( l/«„) (c„= + 2c.^ + 4 c. 2* + + 2’“c„- + . . ) 

X {log (4rt,/c,) - log (rt./rt.) - i log (a.,la.O - [ + «„ (62), 

K’-E'^\ (!,/«„) (1 - - rf -'lei-.. - 2"-'c,“ - ...) 

X {log (4«,/ci) - log («, AO - h log (« 2 /« 3 ) (63) 
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for the calculation of the complete elliptic integrals of complementary 
modulus in terms of the a.o.m. scale (ao = 1, ^ ^0- 

The nome q is defined by the relation log (1/^) from which 

it follows tliat 

log = 2 (D"'* log («»/«n+i) - log (4a, /c,) (64) 

which is more convenient for numerical computation than (60). It may 
be noted that the entire series of calculations may be carried out with 
ordinary logarithms to the base 10, enabling q and powers of q to be 
evaluated directly*. 

The complementary nome q' = is obviously related to q by the 

formula 

log q ■ log q = 

SECTION VII 

LANDEN’S SCALE OP DECREASING AMPLITUDES 

(i) CaUnlatim of the Elliptic Integrals of the First and Second kindsy 
modulus k\ 

If wo denote 

sJ{<X>n COS^ ^n) (65), 

the recurrence formula 

tan (ij/n -- = (cv+i/an+i) tan if/n+i (66) 

may bo written in either of the forms 

sin (2 jAh+i - iA„) - (bn/an) sin cos (2 i/^h+i - ^n) = • . .(67). 

From these we deduce 

4* cos — 2A,j + i, A,| bfiCOSljfn ' *’^^n + i^'n + l/^n+l •••(68), 

which give by differentiation 

# 0/^0 = #l/^l = # 2/^2 = . . . - #n/A„ = (69). 

If we construct the a.g.m. .scale (ao= 1, bo = k') we notice that as n 
increases, xf/n tends rapidly to a finite limit. If we start with 
A^, = A (i/^, k'), wo have, since Cn ultimately tends to zero, 

* The limit derived in (60) does not appear to bo f^enerally known, although an 
equivalent formula is given by Legendre {Fonctions ElUptiqnes, t. i. pp. 94, 101). 
Although formula (60) does not appear to have been explicitly stated by Gauss, 
the scries (61) is found among a number of formulae given in his posthumous 
papers (Gauss, Oes. U'erke, Bd. in. p. 377). Neither (60) nor (61) appears to have 
been noticed by Jacobi or by subsequent writers on elliptic functions. 
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f'f'» 

= A;0 = (l/«n) #n/C 0 Sl/r^ -(1/an) log tan (J7r + il/r„) - -.(TO), 

a useful formula due to Legendre*. 

From (68) we easily obtain the relation 

+ byi^n OOS 

which gives 

Andfn “ ^^n^ldll/n+l = - ^»nCOS 

Writing down the system of equations for w - 0, 1, 2, ... multiplying 
successively by 1, 2, 4, ... 2”, adding and integrating, we obtain the 
result 


f^O f^n+l f^O 

A„^i#,m = J(^o* + 2c,^+4c2H...+2«Oj^ #o/Ao 

- (bo sin i/^o + 26j sin + . . . + 2"6 h sin «/^„) (71). 

Adding to each side the identity 

2"-^^ I bn+i cos = bn +1 Sin i/Tn+i + (2^'*’‘ - 1) bn^i sin 

Jo 

and re-arranging the terms, (71) may be written 

r^o r^n+i ^ , 

I - 2""^^ / (An+i - bn^.1 cos if^n+i) #n+l 

Jo Jo 

f^O 

= j- (Co^ + 2ci^ + . . . + 2^‘c„^) #o/^o + fu+i sin \l/n^i + (bi sin \l/i - />„ sin i/r^) 

Jo 

+ 3 (6j sin i/'s - 6i sin i/-,) + . . . + (2”'^' - 1 ) (6„+, sin i/r^+i - 6„ sin <l>„) 

(72). 

It is not difficult to establish the identity 
6 „h.i sin f„+i -l>n sin >l/„ = 2 c „+2 tan ( 2 vt '„+2 - ^«+,)/cos ( 2\^„+3 - 

...(7.3). 

As n increases indefinitely, «,,+! and tend to the same limit a„, 
while c„+i tends to zero. It follows that the second integral on the left- 
hand side of (72) tends to zero, so that we finally obtain, making use 
of (70), 

E (}li, k') = an sin i/f„ + I (<;„- + 2c, “ + . . . -H 2" c„® + . . . ) 

X (l/a„) log tan (]«• + lf„) + 2 (<j/n, c„) (74), 

2 

where for the sake of brevity we denote 



* Legendre, Fonctions Elhptiquea, t. i, chap. xix. 
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a useful fonnula for the numerical calculation of the second elliptic 
integral*. 

Making use of (27), (28), (48) and (70), we obtain from (74) the 

following interesting expressions 

Z(\l/, k') + l7rF(il/, k')l(fCK') = an sin ipn + 2 (j/^„, c„) ...(76), 

2 

or, in terms of u, by (70) we have 

sin if/n = tanh (a^w), 

and Z (m, k') + anuj K' = tanh { 0 ^ 11 ) + 2 (i//n, o„) 

2 

Jacobi’s relation (47), taken with eciuation (76), gives 

(I/O Z {iu, k) - tan cos (2i/ri - - (in sin - 2 (i^n, c,^ (78). 

2 

The addition theorem 

Z{u)-^ Z {;6) — Z{u-¥ v) - k“ sn {u + v) sn u sn v (79), 

with iu for u and ' 0 -±{K-¥iK') gives two interesting formulae 

iZ{iu + K+ iK\ k) = an(l+ sin iA„) + 2 (lA^, o) 

2 

( 1 /o ^ (m - K- iK\ k) = ( 1 - sin «Aw) “ S (i/^,, , c„) 

2 

If we notice that 




I sin //</'„/ A,. = (l/6„) log («„^.,/A„^.,) - (!//-„) log sec (2i/'„+j - 
Jo 

the result of integrating (72) written in the form 

00 

Z(u, k') + a„ujK' = a„ sin ^n + '2. (2’‘'"*-l) (<>„+i sin i/'j.+i - 6,, sin </'„) 

0 

gives rise to a formula for the numerical calculation of ® (?/, k')^ 

( 81 ), 

in which, if desirable, we may write seci/^n ^osh (</„?/). 


(ii) Calculation of the Jacobian Functions sn {u, X'), cn(;/, k'\ dn (w, k'), 
etc, in terms of the ar(fument u. 

Form the a.o.m. scale (</o-I> lh~k') and determine from 
sin ipn tanh (</„//), n being such that Cn is less than the small quantity 
(letermining the order of accuracy of the calculations. 

* The series (71) was obtained by Legendre {Fonctions Klliptiques, 1. 1. chap, xxi.) 
and left by him in this non.con\ergent foim. As far as the writer is aware, the 
convergent senes (75) is new. 
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From the recurrence formula (66), tan(i/rn-i-tA«) = (<?n/Gn)tan^n, 
compute ^n-i • •• ^ 9 , ^ 1 , and finally terms of which 
sn (w, H) = sin ^o> cn (w, X;') = cos ^o, dn (w, ^r') = = cos (2i/^i — ^o) 

(8Si); 

Z (w, i') may then be calculated from (77) and © (w, K) from (81). 


SECTION VIII 


THE HYPERBOLIC SCALE OF DECREASING AMPLITUDES 

(i) Calculation qf the Elliptic Integrals qf the First and Second kinds, 
modulus k. 

If we write for if/ throughout the formulae (65) to (69), (65) becomes 

An = ^/(an® cosh^ ij/n “ Cn* siuh* if/n) (83), 

while the recurrence formula (66) may be written 

tanh (if/n - ^An+i) = (<?«+i/«n+i) tanh i/fn+i (84), 

or, more conveniently, 

sinh (2lAn+l - «An) = (bjan) sinh l/^n, cosh (2lAn+l - •An) = An/an . . .(85). 
As in Section VII (i), it is readily proved that 

#o/^o = #l/Ai = # 2 /Aa = ... = #n/An = (86). 

If we construct the a.g.m. scale (ao = 1, bo = k') and write 

tan <A = sinh lAo (87), 

it is readily proved that 

Ao = A («A, ^)/cos lA and dipo/^o = <^/A (^, k) (88). 

Thus, writing u = i^(«A> k), we have 

sn (m, k) = sin ^ = tanh lAo, cn (u, k) = cos lA = sech lA©! 
dn (u, i:) = A (^, k) = Ao/cosh lA© = cosh (2iAi - iAo)/co8h iA©j “ 

As w increases, the lA's tend to a limit lAn, while c„ tends to zero. 
Thus the result of integrating (86) gives 

f^n 

u = F(}li,k) = (l/an) j #n/co8h^, = (l/a„) tan-’ (sinh ...(90). 

Jo 

If we notice that 


j*'' \d\l/Q = A (lA, k) #/cos* lA 

= dn® (u, k) du/on^ (u, it) = ^ dn* (iu, if) du, 


K 


2 
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we have from (29) 

f\ = ( 1/0 Z(m, i') + u. 

JQ 

It is easily seen that equation (72) holds good when i\l/n is written 
for xf/n^ Hence if we write 

■S' M, ^ - O/. a. A/, tanh (2»/'a-i^0 

2 *cosh(2i/^s-i/^a) ^cosh (2 i/^4-^s) 


+ 2(2»‘-l)c„^a 


tanh ( 2 >/^n-h 2 - <An+i) 


+ ( 91 ), 


’^cosh (2l/^n+3 - ‘An+a) 
and make use of (47) and (48), we obtain after some reductions 
Z (jA, k) = tanh lAo cosh (2iAi - <Ao) - sinh iA» - S' (lAn, c») . . .(92). 

The relations (27) and (28) enable us to calculate the second elliptic 
integral ii^(«A» ^)- 

(ii) Calculation of the Jacobian Functions sn (m, i), cn (m, k\ dn (^^, i), 
in terms of the argument u. 

Compute the a.o.m. scale ao = l, determine from (90), 

which may be written 

sinh lAn = tan {unU) (93), 

n being such that c» is less than the small quantity determining the 
order of accuracy of the calculations. 

We then calculate the amplitudes lAn+i ••• «A 3 , ^i, from the recur- 
rence formula 

tanh (lAn-l - ’An) = (Cn/^n) taull »An. 

Equation (89) then gives 

sn (u, k) - tanh «Ao, on (w, k) = sech lAo] 
dn {u, k) = cosh (2iAi - iAo)/cosh jAo J 

Equations (91) and (92) then enable us to calculate Z(Uj k). To 
compute ® («/, k\ we note that (92) may also be written 
Z{u, k) = sn (m, k) dn (w, ^)/cn (w, k) - an sinh «An 

- 2 (2**+^ - 1) (^n+i sinh sinh «An). 

0 

Integrating this expression with respect to u and noticing that 

I sinh = (,1/bn) log cosh (2i/',+a - ^„+i), 

Jo 

we obtain the following series for ® (w, ^), 

, r© (u, X:)-| , rcosh i 1 1 1 „ „ r - •A.+O") 

^°4®(o,it)J°^‘"4cosh 

( 95 ). 


( 94 ). 
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SECTION IX 


ON THE NUMERICAL COMPUTATION OF THE 
THIRD ELLIPTIC INTEGRAL 

As defined by Legendre, the third elliptic integral may be written 


n, («, (1 + „ siJt) A (^; 1) 

The parameter n is taken as real, and may take any value, positive 
or negative. 

Writing u = jP(<^, k), n =- sn^ (a, k), Jacobi considers the integral 

V /* “^^sn acnadnasn^M^M v ,, 


‘ sii^ a sir u 


Circular Cases, 


„/ V r “^^sn acnadnasn^M^M 

n (u, a) - / — - — = uZ{a) + Hog -) ' . . .(97). 

^ ' fo 1 “ Arsn^asir w v / - + 

Legendre^s integral is connected with Jacobies by the relation 

Ila (w, ky€l>)-u + Bn a/(cn a dn a) . n (Uy a) (98), 

and in this book is differentiated from it by the suffix s. 

In order to include all real values of n in (96), a is not restricted to 

be real, and the treatment is usually divided into four cases : 

Case 1. n negative between 0 and - 

Case II. « negative between - 1 and - oo } Cam. 

Case III. n negative between and ~ 1) ^ 

Case IV. « positive between 0 and + « | 

Case I. n negativsy between 0 and - Hyperbolic case. 

We write = - X;* sn* (a, A;) = - A;* sin* 0, where 0 <a<K . . .(99). 
Prom (97) and (98) we have 
FT / 7 n sn (a, X;) .Z' (a, ^) 1 

(100). 

^ cn (tt, k) dn (a, k) ® L® (« + Uy k)j ^ ' 

Compute the a.g.m. scale (a© = 1> X;'), and starting with <^o = make 

use of the recurrence formula (17), tan (<f>n+i’-^n) = tan We 
then have sn (a,X;)=sin 0o> cn (a, k) = cos <^o, dn(a, k) = cos 0o/cos(<^i~<^o), 
and JZr(a, X;) =Cisin<^ + <? 2 sin<^a+ ... +c„8in<^n+ ..., while a»a = <j^n/2^ 
We thus have 
n8(n, X:, <^) 

=:u[l+ — — — {ci sin + Ca sin + ••• + Cn sin 

L cos 9o J 

^ ta^^oCOS^^^^) 

^ cos ° L® (a + Uy X;) J ^ ' 
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The calculation of u by the recurrence formula (17), commencing with 
4^0 = <l>i presents no difficulty. The most convenient method of computing 
the last term of (101) is to make use of Jacobies series 

0 (a?) = 1 — cos 2ancc + cos 4an^ “* 2^* cos + (102) 

in which we have written The nome q is most conveniently 

oo 

calculated from (64), log g' = 2 log («n/a«+i) - log (4ai/c0. 

If we write the last term of (101) vanishes, and we obtain for 

the complete elliptic integral of the third kind the expression 

Hj (w, A, Jtt) = ^[1 + tan sec cos (^i - c^o) {ci sin sin <^2 + 

... + CnSinc/>„+ ...}] (103). 

We may also make use of the h 3 rperbolic scale of decreasing amplitudes 
and the formulae of § 8. The final results are not, however, so con- 
venient for computation as (101) and (103). 


Case II. n negative^ between - 1 and - 00 . Hyperbolic case. 


We write 

n = - A® sn® (a + iK\ i) = - l/sn® (a, i) = ~ l/sin® tf, where (^<a<K 

(104). 

Equation (98) becomes 


IIa(w, A:, «^) = 8n*o r — 

^ Jo sn®a-sn®M 

We have 


„ — “ n(M,o+jX') ...(105). 
cnadna ^ ^ ' 


Z(a + iK') = Z(a) + cn a dn o/sn a - ^ 77^7 JT (106) 

and 0 (« + iK') - iBH{z) (107), 

where log B-\ir K IK - ^mzIK, 

H{z) is Jacobi’s Theta-function given by the expansion* 

H{z) = 2g^ (sin sin 3a«« + ^ sin ~ sin + . . .) 

(108). 

Writing successively (a -- v) and (a + u) for z in (107), and a + iK for a 
in (97), we obtain 

n (w, o + iK'^ = uZ (a, £) + 2 / cn a dn a/sn a + J log [H{a. - u)IH(a + tt)], 
and hence firom (105), u<a, 

lls(n. A:, <^) = - Sn a/(cn a dn a) [uZ(a^ ^) + 1 log {H{a~-u)IH{a + tt)}] 

(109). 


* Jacobi, Fundamenta Nova^ §§ 62, 68 ; Oes, Werke^ 1881, Bd. 11 . pp. 229-231 ; 
Cayley, Elliptic FuncHonty 2nd Edition (1895), Chap. vi. 
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For numerical calculation, form the A.G.M. scale («©= 1» RRd 

starting with = make use of the recurrence formula (17), 
tan (<^n+i - <#>n) = (bjan) tan 

We then have ana = <^n/2“, and as in Case I, we derive finally 

<A) 

= -wtan<^o sec ^0 cos (^i -<#>o){ci sin -f Ca sin + ... +c„8in<^H+ •••} 

- J tan <^o sec <^o cos (<t>i - </»o) log [If (a - u)/If (a -i-u)] (110). 

Starting with <f>o = ^ and making use of the same recurrence formula as 
above, the computation of u presents no difficulty. As before, the 
calculation of the last term of (110) is best achieved by computing the 
nome g by (64) and making use of (108). 

If we write ^ = i^r, the last term of (110) vanishes, and we obtain for 
the complete elliptic integral of the third kind the expression 
Hs (n, k, 

= -- A'tan <^osec <#>ocos (<^i~<^o) {ci sin + Ca sin + ... +Cn sin <^»+ ...} 

......(111). 

When w>a, it will be noticed that the integral (105) has a singularity at 
w=a, where the integrand becomes infinite. Atthe same time the logarithm 
of a negative quantity appears in (110), since H (a- u) = - H(u- o). 
It is necessary in this instance to define the path of integration, re- 
garding (105) as an algebraic integral obtained by writing z-snu. Near 
the point z = ma this integral is approximately 

IIj = sn’* a/(cn a dn a) . / a - z^% 

which integrated round a small semi-circle below the real axis gives rise 
to the imaginary terra - sn a/(cn a dn a). If we define the path of 
integration of the algebraic form of the integral as consisting of this 
small semi circle and the remaining portions of the real axis, the above 
imaginary term cancels the term - sn a/(cn a dn a) . | log (— 1) appearing 
on the right-hand side of (110) when we write -£f(a-w) = -i5r(w-a). 
With this meaning attached to the path of integration, we may write 
the incomplete integral in the form, u>a, 

n, (w, i, <^) = - sn a/(cn a . dn a) [uZ (a, i) + 1 log {H {u - a)jH (w + a)}) 

( 112 ). 

Numerical computation proceeds in the manner already described. 

Case III. n negative^ between - and - 1. CircuUx/r case. 

First Method. 

We write 

« = - sn* (ia + JT, *) = - F/dn* (a, K) = - ^/(l - sin* B) ...(113), 

whence sin ^ = sn (a, i'), 0 < a < K\ 
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Writing ia + ^ for a in (98), we find 

IIj (n, k,<l>) = u + i dn (a, ky[k'^ sn (a, k') cn (a, k ')] . n (u^ ia + K) 

(114), 

while (97) gives 

n (u, ia + j^T) = uZ{ia + K) + \ log [0 (ia - w -f jST)/© (ia + w + if)] 

(115). 

Making use of the addition formula (79) for Z{u\ we find after some 
reductions 

iZ (ia + K) = iZ (ia, k) + i^m (a, k')/[cn (a, k') dn (a, k')] ... (1 1 6), 
or, from (47), in terms of a real argument 

iZ (ia + K) = Z(a, K) + i iraKKK*) - sn (a, k^) cn (a, k^An (a, K) 

(117). 

Introducing Jacobi’s Theta-function 0i (z) defined by the relations 

@( 2 ; + = (z) (118), 

01 ( 5 ;) = 1 + 2q cos 2anZ 4 - 2 ^ cos + 2(f cos ^a^z + (119)) 

we may write (114) in either of the forms 

n, («. k, 4 >) = 

n, («, k, 4) - p J 




I'ikK' - Y ^ 

For numerical computation we may make use of two sets of formulae, 
according as we employ circular or hyperbolic functions. 

(i) Use of circular recurrence formulae. 

Form the a.g.m. scale (a© = 1, fro = ^0* Commencing with 4>o = make 
use of the recurrence formula (17), tan (</>„+i - <^„) = (fr»/an) tan<^n, and 
thus derive UnU - <^„/2" = 4>. 

Commencing with ^o = 0, employ the recurrence formula (67), 
8in(2t/r^+i-</r«) = (fr„/an)sinj/'n, giving «na = logtan(j7r + ^i/r„). The 
application of (47) and (76) to both (120) and (121) gives us the formula 

( 122 ) 

for the numerical calculation of the incomplete integral, and the formula 
n, (», *, |>r) = 2 (JST/Ar'*) cos ( 24 i - ^«)/sin 2iA„ . [a, sin + S (^)] 
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The peculiar difficulty arising in the numerical computation of the 
circular cases of the incomplete elliptic integral of the third kind lies in 
the evaluation of the Theta-functions of complex argument. Writing 
= in Jacobi’s expansion (119), the series may be resolved into 
real and imaginary parts*. If we write 

tan Sk ~ sin^2a»t« . sinh 2anO- 2q* sin 4a„^^^inh 

~ 1 + cos 2anU . cosh 2a»a + 2^ cos . cosh 4a»a + . . . 

(124), 

we easily find log [©i - 2 a)/©i + ea)] = 2/^, from which it follows 

that (122) may be written 

Ilg (w, kj </>) = (2lk'^) cos ( 21/^1 - i/'o)/sin 2i/^o {u {a„ sin + S (i^)} - 

• (125). 

Powers of the nome q arising in (124) are most conveniently calculated 
from (64), 

oo 

log = 2 ay'" log (ajan+i) “ log (4«,/cj). 

The use of hyperbolic functions in (124) may be avoided by noting 
that if we write 

Q 2 = ^ + J»A„) = (1 + sin xl/n)l(l - sin i/r,.), 

2 cosh 2a„a = Q® + 2 cosh 4an« - Q* + Q*"*) . 

2 sinh 2a»a = Q* - Q~^) ’ 2 sinh 4a„a = - Q-*) > ' * * ® * * 

from which it follows that (124) may be written 

f Ur - + ((/ - Q -*) sin 4<l> ^ y^(( y - Q' ^) sin 64> + . . . 

^ I + q{Q^-^ Q~^) cos 24> + ^( ^ + cos 44> + q^((/^+ cos 6^ + . . . 

(124 bis). 

(ii) Use of hyperbolic recurrence formulm. 

Write » = - sn^ (w + A) = ~ ^^/(l - tanh* <#»o') (126), 

whence sn (a, k') = tanh <^o', 0 < a < 

From the hyperbolic recurrence formula (37) 
tanh (</»'nfi - </>n ) = (5«/«n) tanh 

we have 

a = <^» 72 ’* = sn (a, k) = tanh <^o', 
cn (a, A') = sech dn (s *0 = sech (</>/ - <^o')- 

* This prooedare is due to Legendre, Fonctioru Elliptiqtust Deuxidme Supplement 
(1828), § iz. p. 146 et $eq. 
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Thus from (44) and (120) we obtain 

tt / 7 1 C08h</>o' 

Ilg (n, k, ft>) - , ,2 , , , . , 


tanh 4>o cosh(<^/ - <^o0 

To calculate u we make use of the recurrence formula (17) 
tan (<l>n+i - 4>n) = (bjan) tan giving anU = <l>nl2^ = 

Equation (124) then becomes 

^ _ 2q sin 2^ sinh 2^' -f- 2(f sin 4^ sinh 4^' + 2g^ sin 6^ sinh 64>' + . . . 


tan^ = : 


1 + 2^ cos24>cosh2^' + 2 ( 2 ^cos 4 ^cosh 4 <^'+ 25 ^®cos 6 <^cosh 6 <^'+ ... 

(128). 

Equation (127) then gives 
v\ f z. A\ ^ cosh (f^Q 

* “ Ftanh<^o'cosh(<^/-<^o') 




and 

n, (», A:, ^ »r) = p -SCnSinh.^* ] 

(129 Us). 

Second Method, 

If we write p = K' -a, (113) takes the form 
w = -~A;2sn®(jr+ziir'--i)8) = -dn®(AA:') = -l+A:'®sin®^, 0<p<K' 

(130). 

We easily find 

n, («, k,i>) = u + P an cn\/9, Jt') 

f- -iri •a\ i -i ® {K iK' - - u)~\ 

X \tuZ{K^rvK -,^)+ J 

dn(A^0 

^ ^ X:'* sn (/8, H) cn (^, A:') 

(131) 

dn(A ^0 

A;'® sn (A A:') cn ()3, Iz) 

r wAr* sn ()&, A:') • fwr-o\ i • i qon 

L" cn (/3, it') dn 03, it') ^ (« - *i3)J “^^^^^’ 

where Hi (z) is Jacobi’s Theta-function defined hy the series 

Hi («) = 2gi (cos a,» -t- gi* cos 3a,» -h g* cos So,* -i- g“ cos 7a,» -»-...) 

(133). 
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Writing ^ f log [H, (u + tp)/H, (134), 

we find 

tan = sin sinh^n P + y^sin 3anU sink San P + y^sin 5anfi sinh ban P-^ • • ♦ 

cosa»wcosha,ji8+g'‘-*cos3an««cosh3a»/S+^cos5a»MCOsh5aH7^+ ••• 

(135). 

The numerical computation may be carried out by two methods. 


Starting with i/^o = ^ in (130), make use of the recurrence formula (67), 
sin - 1 /^„) = (6n/an) sin From (80) we have 

iZ {K + iK' — ip) = (1 - sin i/^n) ~ 2 (j/r), 

while 

sn (A = sin cn ()3, i') = cos , dn (ft A') = cos ( 21 ^^! - i/r J. 
Thus (131) becomes 

n. in, k,4,) = u^ sin + 2 m . . .(136). 

As before, « is most conveniently calculated from the recurrence 
formula tan (<^n+j - <An) = (^n/On) tan <#>„, giving «,« = = *. 

The use of hyperbolic functions in (135) may be avoided by writing 
log O' = <i„/3 = log tan (;i7r + J-i/r„), so that 

2cosha„^=Q'+Q'-M 2 cosh 3a„)8 = Q”’ + O'-’! . . 

2 sinh a„/3 = O' - Q'->J ’ 2 sinh 3a„j8 = O’’ - (/-*J ’ ® 

Thus (135) may be written 

tan sin «!» + !f (Q'» - q-») sin 3<t + / (Q'» - Q'-» )sin5«l>+ - 

(^' + Q'"') cos <*> + (^'» + ^-*) cos 34>' + g* (g»' + ^'-»)co85<t + • • • 

(138). 

Writing m = A'in (136), ^ = Jv, and we have, since iir = a„ir, 

n3(«,^,Jw) = A:[l + 2-‘’||ifc^>{a4l-8in^0-2(^)}]..^ 

(ii) Use of hyperbolic recurrence formulae. 

Write sin ff = sn (ft i') = tanh <^o'> (130) becomes 

w = - A;" sn* (iT + iK' -//?) = -!+ A:'^ tanh* • • .(140). 
Make use of the hyperbolic recurrence formula (37) 
tanh (</)'n+x - ihn) = (ftn/«n) tanh <#>« , 

and we have 

anP = = ^n 72 ^ sn (ft i') = tanh^o', 

cn (ft A;') = sech dn (ft A:') = sech {thi ■“ <Ao')- 

Also -fZ(fft A:) = Ci8inh^' + C28inh<^'+ ... +Cn8inh<^« + ..., 
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SO that (132) gives 

Ha (w, k, 4>) 

1 cosh <t>o 


cosh tanh 

X ui^ sinh <I>q cosh (0/— <l>o) + 2 8inh</>„' + . . .(141). 

As before, we calculate u from the recurrence formula 
tan (<^n+i - <l>n) = (bn/an) tan 
giving = <^>. Equation (135) gives for 

tan ^ + sinh 30' sin 30 + sinh 50' sin 50 + . . . 

cosh O'cos O + (/^cosh 30' cos 30 + g®co8h 50' cos 50 + . . . 

(142). 

Writing u = IC in (141), we have 0'= Jtt, and since 
n, (n, k, ^tt) 

" *'» cosh tanh W f c,8inh<^„'-^8mh<^,'cosh(^,'-<^,')j 

(143). 

Case IV. n positive, between 0 and oo . Circular case. 

First Method, 

Write 71 = ~ ^ sn^ (/a, k) = 1^ sn* (a, i')/cn* (a, k') = ^ tan^ . . (144), 
giving sin^ = 8n(a, A'), 0<a<A^'. 

Thus 

(145). 

®(z) is Jacobi's Theta-function defined by the series 

® ( 5 ;) - 1 - 2^ cos 2anZ 4- 2^ cos 4anZ ~ 2^ cos 6a»« + . . . (146). 

Writing O = Je log [0 (u - ia)/@ (u + /a)] (147), 

we find 

2g sin 2anU sinh 2a„a ~ 2g * sin 4an u^ sinh 4ana + . . 


tanO: 


1 - 2g cos 2a« u cosh 2ana + 2gf* cos 4a» m cosh 4a„a - 


.(148). 


The numerical computation may be carried out by two methods. 

(i) Use of circular recurrence fcyrmulae^ 

Make use of the recurrence formula sin( 2 i^n+i“^«) = (^n/an)sini/r,j, 
commencing with xl/o = 0 in (144). We then have 

sn (a, i') = sin i/^o, cn (a, k') = cos ij/Q, dn (a, k') = cos (2il/i - i/^o)i 
and by (78) 

iZ (ia, k) = an sin + 2 (<A) - sn (a, k') dn (a, X;')/cn (a, k') ...(149). 
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Hence (145) may be written 

(150). 

As before, u is calculated from the recurrence formula 
tan (<^„+, - <^„) = (6n/a„) tan <^n, 

giving a„«=^,/2’*=4>, while a*a = log tan (;!«•+ JV'n)- We have 

2c08h2a„o = Q*+Q-M 2 cosh 4a„a = Q* + I , . 

2sinh2a,.o= ’ 28inh4a„a= 

where Q = tan (Jir + 1 ^„), 

so that (148) becomes 

tan Uf= 5’ (0^ - ^'') sin ^ (Q« - Q -*) sin 44> (^ - (g~«) sin6»-... 
1 ~5'(^+ Q~^)cos 24> + 5 ^ (^+ ^'*)cos 4* - (Q* + Q”*)cos64>— . . . 

„ . (152). 

Writing M = A" in (150), 4' = 0, and 

n, (n, r —1- r - H- «. sin f , + 2 (^)1 

^ ^ ^ 2 COS (2i/^i - i/'o) L tani/fo 2 ^ ^ J 

( 153 ). 


(ii) i/50 0 / hyperbolic recurrence formulae, 

W rite sin ^ = sn (a, i') -= tanli </>o' (1 54), 

so that n--]^ sn* (ta, A:) = F tan* ^ sinh^ 

Starting with given by (154), make use of the recurrence formula 
tanh (<^ n+i - ^n) = (pnM tanh 4>n (155). 

We then have 

and = sn (a, k') = tanh <^o', 

cn (a, k') = sech <l)oi dn (a, k') = sech (<^i' - <^o'), 

00 

and - iZ (ia, A) = X c„ sinh <}>n'. 

Equation (145) may then be written 

D.(%». 56). 

As before, u is calculated from the recurrence formula 
tan - <t>») = (5„/an) tan <^„ 
giving a»M = ^n/2*=#. Equation (148) then gives for ♦ 

2g sinh 24>' si n 2<t - 2g/* sinh 4^' sin 44> + 2</* sinhj6*|^in^4> - 

® ~ 1 — 20 cosh 24>' cos 2 ^+ 2 cosh 4'b' cos 4^— 2 ^'* cosh 6^ 'cos 64' + ... 

(157). 
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Writing u = ^ in (156), = 0, and 

n. (n, k, Jx) - tanh.^.>^coAX^' -<^.0 1 ^ 


Second Method, 


If we write /3 = ^'-a in (145), we have 

w = - A;® sn® (iS7 - i/3) = cn® (i3, A:')/8n* (/3, Ar') = cot® 0 . . .(159), 
giving sn(/3,A:') = 8in^, 0<P<K\ 

We easily find 


Hg (w, At, <^) 

-— ('«»• 

where .£^(5?) is Jacobi's Theta-function defined by 


H(z) = 2q^ (sin sin 1 sin ha^z ~ sin 7an5^ + • • •) 

. (161)- 

If we write 

= |i log [H{u + iP)IH{- u + i/3)] (162), 

we find 

cot -f ^ cos 5aw^ sinh 5anP‘- ■ - » 

~ sin OnU cosh sin coshSani^ + ^sin 5an^^cosh5a»/8~ . . . 

(163). 

The numerical computation can be carried out by two methods. 


(i) Use qf circular recurrence formulae. 

Commencing with i/^o= din (159), make useof therecurrenceformula(67) 
sin (2iAn+i - ’An) = {hM sin fn. 

From (78) we have 

- iZ (i/3, k) = tan cos (2^i - ^o) - »n sin - S (lA) • . • (164), 

a 

while 3n(8, A') = 8in^„, cn(/3,A:') = co8^o, dn (^, i') = cos (2f i - 
Equation (160) then gives 

n,(», «A) 

“ (2.Ai-^.)-«»8in ,A.- 5 (^)} + ♦'] 

( 165 ). 

As before, u is calculated from the recurrence formula 
tan (<^n+i - '^) = (5n/o«) tan 
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giving = ^n/2" = while o»)8 = log tan (Jtt + Thus 

2co8ha,^ = Q' + Q'-‘I 2cosh3a,^ = Q'’+Q">I , . 

2 sinh a„i8= ^ ’ 2 sinh 3o,/9= 

where ^ = tan (Jir + ^ f,). 

Thus (163) becomes 

(Of - ^"Qcos » - (^*- Q -*) co^34> + g* ^~*) cos5»-.. . 

° ~ (Q + Q~‘) sin 4> — g“ (^’ + ^■*) sin 3* + g* (O' “+ O'"') sin 5*- . . . 

(167). 

Writing «= .£" in (163), ♦' = ^ir, and (165) becomes 

n, («, k, |ir) = 2 cfa(2i,-% 7) ““ (1-sin ./'„)-2 (i^-)] 

(168). 

(ii) Use of hyperbolic re(Mrrmce formulae. 

Write sin ff = sn (/8, k") = tanh so that (159) gives 

» = -;fc»8n’(iir'-i^) = cot’<?' = l/sinh“^,', 0<P<K' ...(169). 
Commencing with in (169), use the recurrence formula 

tanh (</>',+i - <t>n) = {bja„) tanh 4>„' (170). 

We then have 

Oifi = <^n72" = S“ (fi> *0 = tanh ‘Ao’. 

cn 08, kf) = sech <l>o, dn (fi, k') = sech (^,' - ^,'), 

oo 

and -iZ(t/3,*) = 2cn8inh<^»'. 

Equation (160) then gives 

n3 (», k, ^) = cosh w [“ i -H *'] ...(171). 

As before, u is calculated from the recurrence formula 
tan («^*+i - <#>«) = (&«/a„) tan 

giving «»« = <^n/2" = *. Equation ( 1 63) then gives for the expression 
, ainh COS ^ — g^ sinh 3^' cos 3^ + (y* sinh 5<I>' cos 5*^ - • ■ • 
cot 4' = ^ cosh 34*' sin 3* + g“ cosh 5*' sin 5* - . . . 

(172). 

Writing w= Ain (163), = and (171) becomes 

n. («, k, i^) = A K + i c>«8inh^;] ...(173). 
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SECTION X 

NOTE ON THE CALCULATION OF THE THIRD ELLIPTIC 
INTEGRAL IN TERMS OF THE COMPLEMENTARY A.G.M. 
SCALE 

When k is very nearly unity, it may in some cases be advantageous 
to make use of the more convergent complementary scale {ai- 1, feo'= ^). 
The problem is to express the various functions appearing in the formula 

k) l 
a + u, k)j 
( 174 ) 

in terms of recurrence formulae involving the complementary a.g.m. 
scale. The four cases discussed in Section IX, depending on the values 
of require separate treatment. The reader will have no difficulty in 
adapting the formulae of Section IX to the complementary scale should 
the occasion arise. We may note, in passing, the formula 

© (a - cn {a + u) © (a k'} 

0(a + u) " cn(a~«4)*©{i(« + *^)> k'} 

which will be found to be useful. 


HaK A:, = w + 


cn 




( 175 ) 



SUMMARY OF FORMULAE 


NUMEEICAL CALCULATION OF ELLIPTIC FUNCTIONS 
AND INTEGRALS 

[Sections containing new formulae are denoted by an asterisk •.] 
First Elliptic Integral. « = F(<^, k) = j* » A =F{\ir, k). 

Second Elliptic Integral. ^) = >/(l - sin’ <b)d<b, E=E{^iT,k). 

Third Elliptic Integral. H, {n, k, <b) = (* 77 — — — a ■ - 

^ (1 +w sin’ </>)./( 1- A® sin’ <#>) 

Jacobian Functions. 

sin ^ = sn {u, k), cos <^ = cn («, k), sin’ ^) = dn (u, k). 


Scale qf Aritkmetico-Geometrical Means. 


a, = ^ («» + *«) 
«j = H«i + ^1) 


bo = k' 
bi = Vao^o 
a^bi^ 


Cq — k 

Cl = i K - io) 
c> ~ i (®i “ bi) 


b. 


1. Landen’s Scale of Increasing Amplitudes 
tan (<^n+i - d>t) = (i«/a«) tan 

F{d>,k)= Lt Ar=H»», 

n-».ao " -J 

(A- A)/^=i(c„’ + 2ci’ + 4Ca>+ ... +2"Cn’+ ...), 

Z A) = E{d>, k) - {EIE) F(<b, k) 

= Cl sin <^i + Cj sin ^ + ... +c«sin<^»+ .... 
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*2. Calculation qf sn i), cn («, A:), dn A) t&rim of argument w. 
Calculate = T'an.u^ n being chosen so that Cn is negligible. Obtain 

^0 from recurrence formula sin(2</>„-i-<^n) = “^sin<^n. 

Then sn (w, A;) = sin <^o> on (w, k) = cos ^o» dn {u^ k) = cos sec {<l>i - <^o)» 


*3. The Hyperbolic Scale of Increasing Amplitudes 

tanh (<^n+i “* <An) = (^n/«n) tanh <#»»• Start with sinh = tan 
F{<h,k')= Lt [lK.<An/2“], 

n-^co 


sinh <^o 


+ ^^V(0.A') 


isr 

— (Ci sinh ^i + C3sinh<;()a+ ... + Cn8inh^n+ •••)* 


*4. Calculation of sn («, A'), cn («, k')^ dn (w, A;') m terms oj argument u. 
(1/i) sn (ew, A) = sn (m, i')/cn (t^, A:'), cn (m, i) = 1/cn A'), 

dn (m, A) = dn {u^ k')lcn (u, K), 

Calculate ^n = n being chosen so that c» is negligible. Obtain 
<An-i ••• <Aa> from recurrence formula 

sinh (2<jin-x - <hn) = (^?nK) sinh 

Then sn (u, K ) = tanh ^o> on (w, k') = sech <;^o, dn {u, k') = sech (<;^i — <j^o)* 

5. Calculation oJ K\ E\ and the nome 

* Hyperbolic recwrrence Jormula tanh Xn+i = \/(^»+ 3 /^»+i) tanh 2x». 
F(</>,k')= Lt [l/an-Xn/^^^], tanh xo = \/>fc' sin 

n-»Qo 

Put<^ = Jir. ^irZ’/jSr=ilog(l/g')=Ltr(^)"log— 1 . (Legendre.) 

L ^nJn-^oo 

log 3 = 2 (I)”-* log ^ - log ^ . (Gauss.) 

1 ®n+l Cl 

_ jE>' = 1 1 (1 _ - 2c,» - ... - 2"-»Cn* - . . .) 

. (log |!-log|-ilog5-. 


6. Landen'a Scale qf Decreasing Amplitudes 
sin (2^»+i - ^») = ^ sin . 
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*E{f,k') = a„sintl/„ + ^(co^ + iCi + ... +2”c„“+ ...) 

X ^ log tan ( Jtt + + S c„), 

2 ^ ' COS (ai/'s - «A2) cos (21/^4 - i/^a) 

+ 2/'2»-lW ton(2.|>„^.a -<^„^.i) 

+ 2(2 

*7. Calculation of sn («, k')^ cn (^^, dn (w, >(;') m <?/’ argument u. 
Calculate from sin if/n = tanh (a»«), n being such that Cn is negligible. 
Obtain \lf n-i ••• ^ 2 , ^0 froiii tan = “ tan ipn* Then 

<^n 

sn (m, A;') = sin if/o, cn (m, k') = cos i/^o, dn (u, k') = cos (2i/^i ~ if/o). 

*S. Hyperbolic Scale of Decreasing Amplitudes 

sinh - W = ~ sinh i/r„. 

Ufi 

P{^i ~ tan~^ (sinh i^„), tan if/ = sinh t/^o , 

jjj 

k)--jy F(il/y k) = tanh \1/q cosh ( 21/^1 - - an sinh if/n - 2' (<An, O- 

jO. 2 


*9. Calculation of sn ^), cn (^^, dn (t^, ^). 

Calculate from sinh if/n = tan w being such that Cn is negligible. 
Obtain •- *^ 2 , ’Ao from tanh (iA„_i - *An) ~ (cn/«n) tanh i^n- Then 
sn (w, k) = tanh i^g, cn {u, k) = sech 1 ^ 0 , dn (m, ^•) = cosh (2iAi - «Ao)/co8h 


*10. Calculation of Third Elliptic Integral — 4 cases, 
sn (a, k) 




Case I . —k'^<n<Q, 
Case III. - 1 < 

Case IV. 0 < < + 00 . 


Case IL — 00 < < - 1. Hyperbolic Cases. 

Method 1. (i) Circular, (ii) Hyperbolic 
recurrence formulae. 

Method 2. (i) Circular, (ii) H)rperbolic 
recurrence formulae. 


Typical Solution. Case HI, Method 1 (i). 

n = - A®/dn^ (a, k') = - h^l{l - sin® ^), sin 6 = sn (a, k') 0 < a < K\ 

1. Form A.G.M. scale {a^^l, b^-k’). 

2. Derive ««« - ^n/2" = ^ from tan (<A»+i “ ^n) = tan <An, commencing 

an 

with <Ao = 
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3. Derive a„o = logtan(J*r + ^^*) from sin (2^„+i - — sin 

(Xn 

commencing with - 0 , 

00 

4. Compute q from log g = 2 log (a«/a«+i) - log (4ai/ci). 

5. Compute ^ from 

tan 'I' ” sin 2a »^ sinh 2ana + 2g* sin 4^anU sinh 4(^a + 2g^+ ... 

~ 1 + 2g cos 2a„^^ cosh 2ana + 2g^ cos 4a„ w cosh 4a„a + 2g® + . . . ’ 

(Legendre.) 

6. Then, 

Da («, ^ -p f 

n, («. /•. Jar) = gin + 2 C„)]. 
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In the preceding sections only those formulae suitable for numerical 
computation have been given. Many interesting and elegant properties 
of the arithmetico-geometrical scales in conjunction with recurrence 
formulae not specially adapted to numerical calculation are given below 
as problems. In some cases the main steps of the proof are briefly in- 
dicated, and unless otherwise stated the final results are new. 
Properties of the a . g . m . scales («©= 1 , h^- 

1. Prove that = <3fn+i + <?n+i> = <*n+i “ ^^n+i • Heuco derive the series 

~ + Cl ^^2 + ^3 + . . . , 

6o = - Cl + C2 + Cs + ... . Gauss, G, W. Bd. iii. p. 376. 

2. Prove that 

\/ Cn+2 a^^ — ~ ^ (>/^n ^n)> ^n+2 — an ^n+2 

and hence Jan= ^€ 2 - 

\/an= ^/c2-^/c4 - *... Gauss, (?. W. Bd. iii. p. 376. 

3. Prove that c^« + 2c\4.i = 2(a» + a„+i)c„+i, and hence derive the 
expansion 

{K - E)IK = (ao + «i) Cl + 4 (aa + a^ 

+ 16(a4 + a5)c5-f ... +22»'(a2n + a2n+i)c2n+i + .... 

4. Prove that 26n-iCn = ic^n~i-2c’*„, and hence derive the series 

= 4 {hoCi + 46iCa + le^aCs + . . . + 2^'^bnCn^.i +...), 
a^ = ^0^ + 2 (poCi + 61 Ca + ... + bnCn-\-\ + . ..). 

5. From the results of Ex. 4, prove that 

EjK-W + 2 (60C1 + 261C3 + 462C3 + ... + 2"“^6n-.iC„ +...), 

EjK = + 2 (hiC2 + S^aCs + . . . + (2**“^ - 1) 6n-iCn +...), 

and hence the identity 

l=an^ + 2 i (2"-^-l)6n-lCn + 

2 0 

6. From the recurrence formula of Ex. 4 prove that 

2 (60C1 + 61C2 + ...) = JCo* - I (Ci* + Ca® + Cs^ + . . .), 

and hence + ^Cq* - ^ (ci* + Ca* + Cg® + . . .), 

= Go® - iCo® - I (Ci^ + Ca^ + C3* + ...). 

Gauss, G, W. Bd. iii. p. 376. 

7. Derive the recurrence formula 

log (an^iM = log (bjbn-i) - 2 log (bn^ilbn) 
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and hence prove that 

IttK'IK = log ( 4 ao/co) - log (ao/«i) - i log {aja^) - . . . 

- (1/2’^) log («n/«n+l) - • • • , 

IttK'IK = log ( 46 i/Co) + I log (62/61) + I log (63/62) + . . . 

+ ( 3 / 2 ”) log {bn,ilbn)+ .... 
Gauss, G. W. Bd. iii. p. 377 . 

[Note that ^irJSr 7 ^= (1/2”) log (ajcn) as n increases indefinitely.] 

8. Derive the recurrence formula 
(l/2”-0 log (ajun-i) = ( 3 / 2”-^0 log («ft/6n) 

+ (l/2”''0 log (ajbn) - (1/2”) log (a„-i/6n-i). 

Hence prove that 

TtKIK = log (16^060/6*0^) + 32 (l/2”)log («n/6n). 

Jacobi, F. N, § 52 , p. 201. 

9 . Show that the a.g.m. scale (ao = 1, h-k') may be calculated from 
the trigonometrical recurrence formula 

sin = tan^ J 

commencing with sin 60 - A, sin 61 ^ (1 - A*')/(l + A')> 

Hence derive the folloiving formulae : 

(i) Ci= JAtani}6o, C2=lJ:tanl6otan^6i, {j 3 =|^itanl 6 otan J-6itan|62> 
Ci‘‘^ = iA*sin6i, C2‘^ = TV^^sin6isin62, C3^ = /YA^sin^isin^2sin68, 

from which, in general, expressions for Cn and Cn may be formed. 

[Note that c„/a„ = sin 6,^, bnjan ~ cos 6,1.] 

(ii) log (2ir/7r) = 2 log (1 + sin tf») = 2 2 log sec 

1 1 

(iii) log {irl{2Kk')} - ^ 2 log sec 6„_i. 

1 

(i v) I irK'iK - - J log ^ log ( 4 /sin 6^) - 2 (^)’‘ log (1 + sin 
[To derive (iv), make use of formula ( 61 ).] 

10. Show that the a.q.m. scale (ao =1, 60 = A;') may be calculated from 
the hyperbolic scale of decreasing amplitudes 

tanh 2xn+i tanh^ X»» or Xn+i = i log cosh 2xn, 
commencing with tanh 2xo = ^, or Xo = i log {(1 + A)/(l - ^)}. 

Hence derive the following formulae : 

(i) Cl - lA-tanhxo, Cs^l^Ajtanhxotanhxi, C3=^i tanh Xo tanh Xi tanh X2> 

c^=\]^ tanh 2xi , - xV^ ‘JXi taoh 2x2, 

tanh 2xi tanh 2x2 tanh 2x3, 
with similar expressions for Cn and c„l 

[Note that 6n/an = sech 2 xn, c„/a» = tanh 2 xn.] 
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(ii) log (2^/ir) = i; log (1 + tanh 2x»). 

(iii) logW(2A'A;')} = 2 2x«- 

(iv) \-aK’IK^ log g" i = log (2 coth Xo) - 2 ( 5 )"log(l + tamhxn+i)- 

1 

[To derive (iv) make use of formula (61).] 

11. Show that the a.o.m. scale (^o - he calculated from 
the hyperbolic scale of increasing amplitudes 

tanh2i/^„ = tanh^i/^„+i, or = | log cosh 2i/^«+i, 
commencing with tanh 2i/^o = or i/^o ^ 1 log {(1 + k')l{l - A')}. 

Hence derive the following formulae : 

(i) C^i+i/Ch = IJa,, = tailh 2ll/n, ejan - sech 2lAn. 

(ii) log (^KItt) = 2 log (1 + sech 2i/r^). 

1 

oo 

(iii) log‘{ir/(2^^')} = 2 log coth 

0 

(iv) = log as n tends to infinite values. 

[To derive (iv) make use of formula (60).] 

12. Writing J {1 + (1 - show that the recurrence formula 

of Example 11 may be written 

2l/r«+l - = log ‘^+ log 

Hence, denoting j8 = (l - + \/k')y derive the expansions 

(i) log q-^ = log (2/P) + 2 (ly^^ log 4+1 . 

(ii) log q ==- log (}P) + 2 ' (}py -^l^(lPy + (1-/3)^= + . . .. 

(iii) q =}jP + 2 (.] py + 15 a py + 1 50 (^py^ + . . .. 

(See also Example 24.) 

13. Derive the following series for tt from the formulae of Example 12, 

7r = log^r"^ — 2.r^- ISa:^ - where a: -I (2^ - 1)/(2^ + 1). 

[Write k = k' = llj2, tlien K~K' and log q"^ -- tt. Note that 

a?-' = 2314086, log a?-' = 3 141602..., 2rc* = *0000070, 
so that two terms of the above series give the value of tt correctly to 
the tenth decimal. The remarkable approximation represented by the 
first term of the series was first pointed out by Legendre, F. t. i.] 

14. Given the nome q, compute the moduli k, k\ 

[The recurrence formula of Example 11 may be written 

4iA» “ 2^/^n+i = ~ log 2 + log (1 - 

= — log 2 — 1 — . . . 
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Note that when n is sufficiently large Vn+i = t log We then 
have 

^ij/n = 2^ log - log 2 - 

Choose n so that the last term is negligible compared to the first two. 
The repeated application of the above recurrence formula enables us 
to calculate successively ••• <^2 and finally iI/q, in terms of which 
k' = tanh k = sech 2i/'o*] 


15. Derive the differential relations 



for all values of n. 


1 J_ 
2“c«^ 


Gauss, G, W, Bd. ni. p. 380. 


16. From the results of Example 15 obtain directly the formula 

[In terms of the a.g.m. scale («o, 5©, ^o) and the complementary scale 
(«o'> Wi <^o)i the formulae of Example 15 give 

But ao - Uqj bo = Co', Co = bo and ao^ = 5©^ + c©®, from which it is easily 
proved that 

Hence we conclude that ~ d log ^ ^ p) ^ 

We also have the identity ^ d ^ (^) ~ which it 

follows that ^ log ^ =- Cn , Cn being a constant independent of 

2 On 

{any bny Cn), (««> ^n, Cn)- In torms of Gauss' notation for the a.g.m., 
we note from Section III that 

an = M{aoy Co) = 2’*if (an, O = 2^anM{\y CnMy 
while an = M{aoy bo)-M{an, bn)* Thus when n becomes indefinitely 

M{ly Cnlan) = ^frlK = ^irllog{4:lk') as k' becomes small. Thus, in the 
limit as n increases indefinitely, Cn = -i^i so that 
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17. Making use of the result of Example 16, prove that each of the 
ratios in Example 15 is given by 

Lt l(anan).d log {an lan)l 

71-^00 

18. Derive Legendre's relation EjK + EjK' - 1 = \'rrl{KK') directly 
from the properties of the a.g.m. scales (ao= ho-K) and {a^- 1, bd-k), 

[Differentiating the identity log (anbn) - log (an-iin-i) = log(a„/5«) and 
making use of the results of Example 15, we have 
d log (Unbn) “ d log (an-ibn.i) = d log (a«/5„) = 2" {cn^jc^^) d log {aolh)^ 
Writing down similar equations for w~-l, n-2y ... 3, 2, 1 in place of 
n and adding, we find 

d log {anbn) - d log {a^bd) = (2Ci^ + 4^2^ + . . . + 2”Cn“) d log (go/^o) (i) 
and 

d log {an bn) - d log {a^bd) = (2ci'* -f 4C2 ^ + . . . + 2^Cn^) cd^^ d log {adjbd) 

(ii)- 

Since -ad=l and bd = Cq, we have 

d log {adbd) - d log {a^bd) = d log {co/bd) ^ ao\'^ d log (ao/^>o)- 

[Ex. 15 with w = 0.] 

When n becomes large, an -bn, so that 

d log {anbn)-d\og {an bn) = 2d log {ajan) = - (4/7r) anUnCf^ d log (ao/W- 

[Ex. 17.] 

Note also the identity d log (oo/^o) + <^d~^ d log {adibd) = 0. [Ex. 16.] 
Subtracting (ii) from (i) and dividing out by 

Co~* d log {a^lbd) = - cd-^ d log {ad tbd) 

we find 

(4/^)anan - Go® 

= - (2ci" + 4^3^ + . . . + + 4 C 3 '' + . . . + 2’‘Cn'" + . . . ) 

(iii)* 

But we have {K - E)IK = \ {cq + 2ci^ + . . . + 2’‘c»* +...),«» = «o = 1 

with similar expressions for the complementary functions in terms of 
accented S 3 anbols. Thus (iii) becomes 

irl{KK') - 1 = - 2 (1 - EIK) + Co** - 2 (1 - EIK') + Co'** 
or, noting that + cd^ - 1, we derive Legendre's relation in the familiar 
form 

EjK^ ElK^ - 1 = \’^l{KKy 

The equivalent of the series (iii) is given by Gauss (G. W, Bd. in. 
p. 380), but it does not seem to have been hitherto pointed out that 
Legendre's relation follows immediately from it.] 



40 


APPENDIX 


19. Verify the following identities : 

(1 + 20 ^ + 20 ^ + 20 ^+ ...)(1 — 2 .r + 2 ^ — 2 ^+ ...) 

= (l-2^- + 2^-2^«+...)", 

(1 + 2^ + 2^ + 2^*+ ...)^ + (l~2^ + 2^^-2.'r®+ ...)^ 

= 2 (1 + 2 ^ + 2 ^ + + . . . )% 

(1 + 2^ + 2^-^ 2^ + -2^ 4- 2^-2^+ ...y 

= (2oi + 20^ + 2x^' + .../. 
Gauss, NachlasSy G. W. Bd. iii. p. 447. 

20. From the results of Example 19, prove that if we form the 
A.G.M. scale commencing with 

ao = h(l +2o + 2o^ 4- . . 6o - ^ (1 — 2^2? 4- 2o* — • • 0^ 

Co = A (2.2?^ + 2x* 4- 2o^^ 4- . 

the wth term is given by 

an = /r(l 4- 20n 4* 20n‘^ + • • •)^ = A (l - 2Xn 4- 2Xn^ ~ • .)S 

Cn = A (20’n^ 4- 2Xn^ 4- 2Xn''*' 4- . . .)2, 

where log Xn - 2” log x and A = Micb^y AJ. 

Gauss, NachlasSy G. W. Bd. iii. p. 448. 

21. Prove that x as defined in Example 19 is identical with the 
nome q when (go-I, Iq-K), 

[When n becomes large, 

- ^ log <7 = Jtt /T'/A" = Lt “log(a«^n) 

n-^oo ^ 

"" Lt i log [( I - 2a?„ + . . . fl^x} ( I + 2a?„“ + • • • )T 

n'^oo a 

= - Lt ~loga:„ 

n-^oo ^ 

= -ilOgXy 

SO that x^q.^ 

22. From the results of Examples 20 and 21, derive the following 
expressions : 

Jp~ 1 -2<74-2(/^-2^*4- ... . , _ 2(/i4-2<7® 4-2//^ + ... 

^ ' “ r+2</ 4- 2(fT2q ^+\ . . ’ ^ ■ 1 + 2q^^~+ 2^4-... ’ 

K-\Tr{\ + 2q + 2q* 4 - ...)^ Kk' \Tr{\-2q + 2q^- 
[From Examjde 20 we obviously have = 1, Aq = k\ Cq = A, A = 

23. If we denote ^ = (1- ^/k')/(l + Jk'), derive the following ex- 
pressions : 

iP=^(q + (f + q^+ •••^(1 + 2^^ + 2</^® 4- ...), 

Ar= 27r (1 4- 2q* + 2q^^ + 4- JkJ = Sir(q + q^ + q^+ . ..)V(l - JkJ. 
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24. From the preceding examples derive the well-known expansion 

q = iP + 2 ( + 15 (iPy +150 + . . . etc. (See also Example 12. ) 

25. If we write 

1 + 2.rn + 2^n‘+ ..., /?n 1 - 22 ’h + 

yn = + 2a^J + + ..., 

where logar„ = 2" logg, prove that the recurrence formula 
tan = (a„+j/^„+,) tan 2i/'„ 

leads to the result 

Lt (iA„/2") = (^irllC) k), 

n-»»oo 

where tan j/^o = 'Jk' tan <#>, and K is given by the (/-series of Example 22. 

26. In terms of Landen's scale of increasing amplitud(»s (Section IV) 
tan - <i^n) = (Kl^n) tail <^n, derive the following expansions : 

dn (Uy A) = A (<^o) k) = an + Ci cos + 6*2 cos <^2 + • • • + Cn cos + . . , 

oo 

sn umu = ^ 2 2’‘“^c\_i sin cos (2 </>h-i - <t>n)- 
0 

27. When (c»+i/(«n+i) is sufficiently small, prove that <^n+i of Example 
26 may be computed from the expansion 

- 2<f>n-(Cn^ilan+i) sill 2<l>n 

+ i (CnJan^lf sin 4:<l>n - i (c„+iK+l)* sin 6c^n “ • • • . 

28. Prove that the recurrence formula of Example 26 is equivalent 
to A„ = J (A\_i + b^n)l^n-i, and hence show that 

log [® (0)] = i: -7)^1 log (a„/A„). 

0 ^ 

[Jacobi, Fundamenta Nova (1881), § 52, p. 204.] 

29. Prove that the recurrence formula an tan ff>n+i = ^n+i tan tfin^ where 
= \l{€in COS* + bn siu* <^n), may also be written in the forms 

sin <f>n+i __ tan 4>n + l _ ^n + l _ '^bji 

sin (bn a„ + A,/ tan<^n s/an s/an + ^n' n/^h + ^h* 

30. Prove that as a result of the recurrence formula of Example 29, 
d(bnl^n - d(bn+il^n+i- Hcuce derive the following formulae : 

F{(l>, k) = Lt (^n/a^t), where = ^0, sin (bo = sn (w, k), 

Z {(b, k) = 2oi cos (bo sin (bi + 4^2 cos sin <^2 + • • • + 2”c„ coh ( bn - 1 sin<^,i + . . ., 
log [® (a, A)/® (0, i)] = 2 2" log [2a„/(a„ + A„)]. 
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[The recurrence formula which lies at the basis of the results obtained 
in Examples 29 and 30 was originally given by Gauss in 1818 in his 
paper “ Determinatio Attractionis, etc.” G, W, Bd. m. (1866), p. 353, 
and was further discussed by Jacobi, “Zur Numerischen Berechnung 
der Elliptischen Functionen,” G, W. Bd. i. (1881), p. 351. It is easily 
seen that these formulae are not adapted to rapid numerical calculation.] 


31. If sin <^o = sn {u, k) prove that 

A, A,.. .An- 


tan a^u = - 


tan<i!>o 


a^aia ^ ... 

[The result follows from Gauss' recurrence formula of Example 29 ; 
see also Jacobi, F. AT., §36, G. W, Bd. i. p. 154.] 


32. Prove that Landen's recurrence formula 

tan (l/^n ”■ ’An+l) ~ (^n+l/^n+l) tan */^n+l 
may be written tan 2xn+i = tan Xn if we write 

tan il/n = Jianjc^) tan x«. 


33. If we make use of Landen’s scale of decreasing amplitudes 
sin (2iAn+i - ^n) = (&«/»«) sin 

commencing with prove that K= {^Kjir) log tan (Jtt + Ji/^n) when 

n increases indefinitely. Hence show that the nome q may be computed 

from the formula q = cot^ (Jtt + \\f;^ = (1 — sin <An)/(l + sin j/^„). Show, 

furthermore, that Legendre's relation leads to the identity 

an (1 - sin ^n) = ( WA') {1 - tanh (jrKl2K)} = 5 («An, c,*), 

2 

where 2 Cn) is the series of equation (75). 

2 


34. Prove that if the a.q.m. scale (^o* ^o) is computed, the circum- 
ference I of an ellipse of semi-axes (^o, ^o) is given by 

I = (27r/an) {a^ - 2c^ - Acs* - ~ . . .). 

35. If we denote by Vi and ra the greatest and least distances re- 
spectively of a point on one circle to the circumference of a coaxial one, 
show that in terms of the a.g.m. scale (^0 = W = the coefficient of 
mutual induction is given by 

M- {27r*/a„} [ci + 2^2^ + 4c8* -I- . . . + +•••], 

or, when the circles are very close together, in terms of the comple- 
mentary scale (do' = n + ^ai = ri - ra), by the formula 

[King, L. V., Proc, Boy, Soc, A, Vol. c. 1921, p. 63 et seq. Several 
numerical illustrations are given in this paper.] 
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